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ABSTRACT 
The  pressure  and  phase  distribution  of  sound  in  a  fast  fluid 
medium  underlying  a  tapered  fluid  medium  was  modeled  using  a  Green's 
function  approach.   The  model  predicted  well  defined  beams  in  the 
bottom,  with  apparent  interference  effects  strongly  evident  for  cases 
at  a  distance  from  the  wedge  interface.   The  effects  of  attenuation  on 
the  patterns  were  studied.   A  simple  expression  for  determining  the 
beam  depression  angle  as  a  function  the  wedge  angle,  the  sound  speed 
ratio,  and  the  density  ratio  was  derived.   Comparisons  with  laboratory 
measurements  were  made. 
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I.   INTRODUCTION 

The  problem  of  radiation  behavior  in  a  wedge-shaped  medium 
overlying  a  bottom  possesing  a  greater  speed  of  sound  has  been 
investigated  both  optically  and  acoustically.   In  1971,  Tien  and 
Martin  (Ref.  1)  deposited  a  thin,  tapered,  dielectric  film  on  a 
substrate  having  a  higher  refractive  index,  and  observed  the  behavior 
of  a  laser  beam  coupled  into  the  film.   Perfect  reflection  at  the 
film- substrate  interface  was  observed  until  the  changing  angle  of 
incidence  decreased  below  the  critical  angle  for  continued  mode 
propagation.   For  angles  of  incidence  below  critical,  the  light  was 
converted  into  radiation  modes  in  the  substrate. 

A  concurrent  analysis  of  the  acoustic  problem  was  done  by  Kuznetzov 
(Ref.  2).   Kuznetzov  developed  a  theory,  based  on  normal-modes,  that 
predicted  the  following  for  sound  propagating  within  a  wedge-shaped 
medium  overlying  a  half -space  with  a  higher  sound  speed: 

1.  As  the  sound  approached  the  vertex  of  the  wedge,  it  would 
continue  to  be  totally  reflected  until  the  angle  of  incidence  decreased 
to  the  limiting  angle  of  total  reflection. 

2.  Sound  incident  at  less  than  the  limiting  angle  of  total 
reflection  would  be  totally  refracted  into  the  underlying  half-space. 
This  refraction  occurs  within  the  region  from  the  wedge  apex  to  the 
point  where  the  limiting  angle  was  attained. 


3.      Within   the  half-space,    the  acoustic  energy  would  be  propagated 
as  a   well   defined  beam  with   the  maximum  pressure  amplitude  occurring  at 
an  angle  of  depression   from  the  wedge  interface  within   the  range   6  to  23, 
where   8  is   the  wedge  angle. 

Kuznetzov  also  performed  a   series  of  experiments   that   supported  his 
theory . 

Subsequent  optical   experiments  by   Tien,   Smolinsky ,   and  Martin    (Ref .    3) 
demonstrated   the  formation  of  a   well-defined  beam  within   the  substrate. 
The  accompanying  theoretical   discussion  predicted  conflicting  behavior. 
Ray-optics  predicted  refraction  into   the  substrate  after  the  cutoff  dis- 
tance   (the  point  where  refraction   begins)    and  wave   theory  predicted  com- 
plete refraction  into   the   substrate  before   the  cutoff  point.      Sigelman, 
et.    al.     (Ref.    4)    both  predicted  and  observed  a   well   columnated  beam  in  a 
water-aluminum  wedge  system,    but   the  observed  beam  was  much  broader   than 
predicted.      Also,    the  pressure  maxima   for   the  columnated  beam  occurred  at 
approximately   11°   for  a   wedge  angle  of  1.3° ,    well   outside   the  limit  pre- 
dicted be  Kuznetzov. 

The  primary  objective  of  this   theoretical   development   is   to  produce 
a  simple,    fast   computer  model    to  predict   the  pressure  and  phase  distri- 
bution of  the  columnated  sound  beam  in   the  fast  substrate .      The  pressure 
and  phase  distribution  along  the  wedge  interface  calculated  by  Kawamura 
(Ref.    5)    will   be  used  as  input   to   the  substrate  model,   and  will   specifi- 
cally be   used  as   the  boundary  conditions  for  a  Green's  function  analysis . 
To  remove  explicit  frequency  dependence,   all   relevant   distances  will   be 
normalized  to   the  distance  measured  from  the  apex  along  the  wedge  inter- 
face at  which  the  critical   reflection  angle  is   first  exceeded;    this    dump 


distance  is   referred  to  as   the  dump  distance  X.      Integration  over   the 
wedge  interface  will   be  truncated  at  distances  sufficient   to  encompass 
only  selected  modes  of  interest.      Results  will   be  analyzed  for   trends 
and  compared  with  existing  measured  data. 


II.   THEORY 

Pressure  and  phase  distributions  along  the  bottom  of  a  wedge- 
shaped  medium  are  obtained  using  an  adaptation  of  the  method  of  images 
as  prescribed  in  Ref.  5.   The  source  may  be  set  at  either  finite  or 
infinite  distance,  but  must  be  maintained  sufficiently  far  from  the 
wedge  boundaries  to  ensure  that    impinging  waves  are  essentially  locally 
planar.   The  following  theoretical  discussion  uses  the  above  pressure 
and  phase  distribution  as  initial  conditions,  and  calculates  the 
pressure  and  phase  distribution  in  the  fast  bottom  by  utilizing  the 
Green's  function  (Ref.  6). 

Figure  II-l  illustrates  the  geometry  of  the  problem.   The  horizontal 
x-axis  and  the  vertical  z-axis  have  their  point  of  origin  at  the  wedge 
apex,  with  the  positive  x-axis  forming  the  interface  between  the  wedge 
and  fast  bottom.   The  wedge  angle  is  designated  S  and  the  angle 
specifying  the  sound  source  location  is  6 .   X  is  the  first  mode 
dumping  distance,  and  is  a  distance  determined  as  follows: 


9   =  arccos  (  C2  /  Ci  )  (Il-la) 
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D   =  Ci  /  {  4«f«sin(0  )  }  (Il-lb) 
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where 

Ci  =  sound  speed  in  the  wedge 

C2  =  sound  speed  in  the  bottom 

f  =  frequency 


and  D  =  lowest  mode  cut-off  depth. 


To  make  the  problem  frequency  independent,  all  distances  were 
normalized  to  X  by  defining 

I   =   x  /  X  (II-2a) 

V   -  x'  /  X  (TI-2b) 

and  n  =  z  /  x  Cll-2c) 

The  variable  x'  is  identical  to  x,  and  g'  is  identical  to  £,   The 
prime  designation  is  required  for  dummy  variables  of  integration  in 
the  Green's  function  development  to  follow. 

Pressure  along  the  wedge  bottom  is  then  given  by 

p(C,n=0,t)  =  p(0  exp  {iCoat  +  L(E)]>  ClI-3) 

where  P ( ^ J  is  the  pressure  amplitude  along  the  wedge  bottom  and  LC£5 
is  the  phase  along  the  wedge  bottom.   Pressure  in  the  bottom,  away 
from  the  boundary,  is  then  given  by  Ref ,  7, 

p(r,t)  =  exp[iu>t]  7-  (  P(^)  |fr     X  dV  (n-4) 

4lT  I  3Z  'z'=0 

=  exp  [iwt]  R(r) 
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where  G  is  the  appropriate  Green's  function.  The  z'  coordinate  is 
identical  to  the  z  coordinate  and  r  =  (x2  +  z2 )  .  The  appropriate 
Green's  function,  for  a  line  source  in  cylindrical  geometry  is 

G(r,r')  =  1ttH(2)  (k2|r  -  r- I )  (II-5) 

o 

where  k2  is  the  propagation  vector  in  the  bottom  and  r  =  [ (x" ) 2  +  (z" ) 2] 

Differentiation  of  the  Green's  function  yields  (Ref.  7) 

^  ,    k2  (z  -  z') 

t-t-  =  1ttH<2;  (k2  r  -  r*  )  (II-6) 

I  r  -  r  I 

where  k2 1  r  -  r   >>  2ir  and  H0  ,  H,  are  zero  and  first  order  Bessel 
functions  of  the  third  kind  (Hankel  functions)  respectively.   Using 
equations  (II-2a,b,c)  at  z"  =  0  the  expression  [r  -  t"\      may  be 
rewritten  as 


r  -  T~\    m   x 


[(?  -  D2  +  r)2]h  CII-7) 


3G 
Substituting  (II-7)  into  equation  (II-6) ,   and  evaluating  - — —  at  z'  =  0 

3  Z  i 


yields 

3G_ 

3z' 


vJzfcA     *xp{f-)  exp{-ik2x[(g-r)2   +  n2]ia>k2nx 

«--o "  x^[(5-r)2  +  n2]3/H 

Combining  equations  (II-8)  and  (11-4)  yields  a  time  independent 

-*■ 
expression  for  amplitude  and  phase,  R(r) 

1/2 


(II-8) 


R(r)    =^exp( 


^m 


rP(V)eXp{i[L(V)-k7x[&-V)2+^2]1/2l}dV 

I  [(K-V)WY/4 

Equation    (II-9)    is   the  basic  equation   used  for   the  calculation  of 
phase  and  amplitude  in   the  bottom. 
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III.   DEVELOPMENT  OF  WDGBTM 

Program  WDGBTM  was  designed  to  perform  four  main  functions. 
First,  the  pressure  and  phase  along  the  interface  of  a  wedge  and 
underlying  fast  media  may  be  calculated  in  two  ways,  for  both 
infinitely  and  finitely  distant  sound  sources  in  accordance  with 
Ref.  5.   Modifications  were  made  to  the  code  of  Ref.  5  to  facilitate 
subsequent  in-bottom  calculations,  and  to  enhance  user  flexibility. 
These  changes  did  not  affect  results.   Then  using  the  first  code  as 
input,  additional  code  enables  both  the  computation  of  pressure  and 
phase  along  a  line  perpendicular  to  the  wedge  bottom,  and  as  a  function 
of  radius  and  angle  measured  from  the  first  mode  dumping  point,  X. 

Modes  of  program  execution  may  be  chosen  by  literal  input  directives 
combined  with  appropriate  numerical  input.   Appendix  A.  illustrates  an 
input  set  that  utilizes  all  four  options.   Output  consists  of  both 
printed  listings  and  graphics. 

Program  versions  exist  for  execution  on  both  a  CDC  6500  and  an  IBM 
360.   The  CDC  6500,  with  a  60  bit  word,  does  not  require  double 
precision  for  accuracy  equivalent  to  the  IBM  360  output.   The  CDC 
version  of  WDGBTM  has  simpler  varian  control  routines,  and  separate 
output  files  for  each  pressure/phase  variation.   This  enables  optional 
suppression  of  listings.   Execution  time  on  the  CDC  6500  is  roughly  ten 
times  faster  than  on  the  IBM  360.   For  average  input  sets,  times  are  150 
and  1500  cp  seconds  respectively. 
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A.   EQUATION  PREPARATION 
1.   Depth  Domain 

Integrations  along  lines  perpendicular  to  the  wedge  bottom 
were  designated  depth  domain  calculations.   The  program  was  designed  to 
allow  specification  of  the  field  point,  i.e.,  specification  of  where 
along  the  media  interface  line  the  perpendicular  probe  line  could  be 
dropped.   The  field  point  is  not  confined  to  lines  directly  beneath  the 
wedge,  but  may  extend  to  points  beyond  the  wedge  apex.   These  points 
would  be  values  of  negative  x   in  Figure  II- 1.   The  beginning  and  ending 
points  of  the  section  of  the  perpendicular  to  be  examined,  and  the 
number  of  field  points  used  in  the  interval  were  also  to  be  specified. 

Equation  (II-9)  was  split  into  real  and  imaginary  components  before 
coding.   Integration  order  selected  was  over  the  entire  £ ^domain  for 
each  H  point. 

2.   Angle  Domain 

Examination  of  Figure II-l   verifies  the  following  geometric 
relationships . 


e  =  a  cosQ7  (Ill-la) 

k 


and 


n  =  (o2  -  z2)1/2  (Ill-lb) 


where  £  is  the  distance  from  X  to  the  %'   coordinate  of  the  field  point 
of  interest  and  C  is  the  radial  distance  from  X  to  the  field  point. 
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Substitution  of  equations  (Ill-la)  and  (Ill-lb)  into  equation  (II-9) 

yields 

2      2    1/7 


p,t,    -    (02-Z2)1/2  /  iTT  \    I2k2x\1/2 

R(r)    = exp     —       — £-  (III-2) 


(*P  (V  )  exp{i[L(V  )  -k2X[(Z,-V  )2 +g2  -z2]1/2]}dV 

Equation  (III-2)  is  the  basis  for  angle  domain  calculations.   The 
specification  of  a  originating  at  X  was  chosen  to  conform  to  existing 
experimental  data  (Reference  8) .   The  length  of  the  radius,  the 
beginning  and  ending  of  the  angular  sweep,  and  the  angular  increment 
comprise  the  input  specifications.   As  in  the  depth  domain  option, 
equation  (II 1-2)  was  split  into  real  and  imaginary  components  before 
coding. 

3 .  Design  Restrictions 

Theoretical  development  of  program  WDGBTM  assumes  the  sound 
source  is  sufficiently  distant  to  ensure  essentially  locally  planar 
waves  at  the  interface  and  the  in-bottom  measurements  satisfy  the 
condition  k  |r-r|>>2iT. 

4 .  Program  Results 

An  initial  set  of  WDGBTM  runs  were  made  in  the  depth  domain  for 
a  perpendicular  directly  down  from  the  wedge  apex.   The  source  was 
infinitely  distant  and  set  at  an  angle  of  6/2.  This  limited  modes  within 
the  wedge  to  odd  numbered  (symmetric)  modes  only.   Integration  over  the 
wedge  bottom  was  truncated  at  5  equal  four,  so  that  only  modes  one  and 
three  were  considered  as  contributing  significantly  to  the  beam  pattern 
in  the  bottom.  Beam  patterns  in   the  bottom  were  calculated  over   the  range 
n  equal   zero   to  1 .  3X. 

15 


The  initial   set  was  comprised  of  seventy-five  runs.      In  addition 
to  the  above  criteria,    the  runs  were  generated  by  all  possible  combi- 
nations of  the  following  specifications.      The  ration   p  /p  was   0.95, 
0.90,    0.80,    0.70  or  0.50.      The  sound  speed  ratio,   C  /C     was  set  at   0.95, 
0.90,   or   0.85.      The  wedge  angle   6  was   varied   through  10,    6.92,    5,    2   and 
1  degree  (s) . 

From  the  results  of  the  initial   data  set,    the  depression  angle  and 
beam  width  was   calculated  for  the  beam  pattern  in   the  bottom.      Results 
are  summarized  in  Appendix  B.      These  results  were   then  used  to  find  a 
simple  analytical   formula   to  predict   the  depression  angle  as  discussed 
in  Section  IV. 

Other  runs  were  made  for  perpendiculars  at  x  =  -9X  from  the  wedge 
apex.      These  runs  displayed  apparent  strong  interference  phenomona  and 
are  discussed  in  detail   later. 

Angle  domain  executions  were  made  to  compare  with  experimental   data. 
Results  are  discussed  in  Section   VII. 

B.       SENSITIVITY   TO   INTEGRATION   STEP   SIZE 

It  was  desirable   to  keep  both  the  integration  step  size  and  the  field 
point   increment   size  as   large  as  possible  in  order  to  minimize   total 
computer  resources  required.      Increment  of  0.01    times   the  dump  distance 
(X)   proved  to  be  optimum  for  both  integration  increments  along  the 
wedge  interface  and  field  point  steps  in   the  bottom. 

To  illustrate   the  adequacy  of  0.01X  increment,   a  series  of  runs  were 
made  in  accordance  with   Table  III-l.      An  overlay  of  wedge  interface 
pressure  amplitude  for     Ax  =  0.01X  and   Ax  =  0.001X  demonstrated  no  dis- 
tinguishable difference.      Likewise,    depth  domain  plots    (a,    b,    c,   d  of 
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Table  (III  -  1) 
Sensitivity  to  Integration  Step  Size  (Ax)  and  Field  point  interval 
size  (An) 

Ax  =  0.001  Ax  =  0.0001 


bvf=0.01      t\r\=0.001 


t\r\=o.oi    An=o.0Oi 


(a) 


(b) 


(c] 


(d) 


Run 

a 
b 

c  " 
d 


H  amp  max 

0.240 
0.244 
0.240 
0.243 


e 


D 

13.50 
13.71 
13.50 
13.66 
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Table  III-l)    could  be  identically  superimposed.      Use  of  an  integration 
step  size  of  O.OOlX  resulted  in  a   calculated  depression  angle  within 
one  percent  of  the  depression  angle  for  the  O.OlX  run.      Utilization  of 
a  step  size  of  greater   than  O.OlX  is  not  recommended.      Within   this 
limitation   the  integrations  can  be  carried  out  over  10X  along  the  inter- 
face for  a  maximum  field  interval  of  10X  in   the  bottom.      The  10X  re- 
striction is   imposed  by  program  array  size.      (For  the   o'-s  specified,    0.5 
degrees  in  the  angle  domain  was  slightly  finer   than   the  O.OlX  field 
point  step  size  in   the  depth  domain  and  was   the  field  point  increment  of 
choice) . 


IV.   BEAM  DEPRESSION  ANGLE  ANALYSIS 

A  simple  relationship  was  sought  between  the  physical  parameters 
associated  with  the  wedge  system  and  the  resultant  beam. 
Assume 

where  9  is  the  beam  depression  angle  and  3  is  the  wedge  angle.   For 

.L/C 

each  p . /p .   and  C./C   ,  9^ vs  B  was  graphed  on  log-log  paper,  and  the 
slope  (  n)    was  measured.   Results  are  listed  in  Table  1  of  Appendix  C, 
n     equals  0.329. 

Then  9   was  found  as  a  function  of  C../C  .   The  critical  angle,   9  , 
is  an  explicit  function  of  C./C 

9  =  arccos    (C  /C  )  ,9C  in  radians  (IV-2) 

It  was  then  assumed 


v-^2iV 


For  any  constant  S ,  M   was  derived  as  the  slope  of  a  log-log  graph  of 
9   vs  9  .   Graphs  were  analysed  for  6  equal  to  10,  6.92,  5,  2  and  1 
degree,  and  P-,/P2  ratios  of  0.95,  0.90,  0.85,  0.70  and  0.50.   Results 
are  summarized  in  Table  (2)  of  Appendix  C,   M  =   0.772. 
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The  same  process  was  repeated  for  ascertaining  dependence  on 
Pl//P2  *   Lo9-l°g  graphs  were  prepared  by  holding  8  and  9  constant 
using  the  assumed  relation 


9dc  -  e0-329  bc0-772  ^/P2>q  dv-4, 


Results  are  summarized  in  Table  3  of  Appendix  C  q  -   -0.254. 

The  proposed  formula  is  now  dependent  on  all  physical  properties, 
assuming  the  source  angle  is  maintained  at  one-half  3.  To  remove  the 
proportionality  sign,  write 

n0.329   n    0.772    ,       .        -0.254  ,,   _, 

e^,  -  «  6         ec         rP2/P2;  wv-sj 

and  solve  for  K  .      From  Table  4,  Appendix  C,  K=   17. 22/radian . 

Therefore 


9dc  -  17.22   S°-J29   a/"772    (P^r0-254  (IV-6, 

where  0  is  in  degrees  aid   ®cin   radians. 

9   calculated  from  equation  (IV-6)  and  the  difference  between  9__ 

and  9  from  program  WDGBTM  (labeled  AQ  and  9   respectively)  are  also 
D  D  DW 

tabulated  on  Table  4,  Appendix  C.   A9  varies  from  a  minimum  of  0.00  to 
a  maximum  of  1.25,  with  the  average  of  the  absolute  values  of  A9  equal 
to  0.36.   Errors  ranged  from  zero  percent  to  5.0  percent.  Thus    9   gives 
a  reasonable  approximation  of  the  beam  depression  angle. 
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A  similar  formula 


ea  -  27.2  S2/3  8C0-773  ^/p2r0-27  ^-/; 


was  developed  independently  by  A.  Coppens  (Reference  6)   The  minimum 

deviation  of  9  ,  from  9  is  0.01  and  the  maximum  deviation  is  1.56. 
DA  DW 

The  average  of  the  absolute  value  of  the  deviation  is  0.36.   Both 
formulas  are  best  for  P-,/P9   close  to  1.0,  and  less  accurate  for 
p,/p2   equal  to  0.5. 
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V.   THE  EFFECT  OF  ATTENUATION 

As  previewed  in  Section  III-A.4  the  pressure  distribution  as  a 
function  of  n  at  x     =  -9.0X  demonstrated  apparent  interference  effects 
(see  figure  V-l) .   As  integration  over  the  wedge  interface  to  4X 
encompasses  dumping  by  the  first  and  third  mode,  and  the  third  mode 
dumps  about  the  point  3X  vice  X  for  the  first  mode,  it  is  logical  to 
expect  attenuation  effects  to  reduce  the  contribution  of  the  more 
distant  mode.   The  initial  Green's  function  approach,  equation  (II-9) 
did  not  consider  attenuation  of  the  acoustic  wave  in  the  bottom 
medium.   An  attenuation  term 

exp  {-a[r5-rj2  +  n2]i/2}  (V-l) 

was  incorporated  into  equation  (II-9)  to  yield 

1/2 

R(h    -n/4exp(^)(^)  (V-2) 

/fprcOexp{-a[^-r;2^2]I/2}exp{i[Lrr;-^2^[rc-r;2^n21i/2]}dr 

An  appropriate  code  modification  to  the  depth  domain  routine  was  made. 
The  attenuation  coefficient  a  was  an  input  variable  to  the  program. 

Two  sets  of  comparison  runs  were  made.   The  first  set  consisted  of 
an  unattenuated  control  (01=  0)  and  an  attenuated  run  with  a  equal  1.0, 
with  the  pressure  amplitude  measured  down  from  the  wedge  apex.   The 
second  set  consisted  of  an  unattenuated  control  plus  runs  with  a  set  to 
1.0,  0.5,  0.1,  and  0.01  and  the  pressure  measured  at  X   =  -9.0  from  the 
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apex.   All  other  input  variables  for  both  sets  were  held  constant  and 
are  specified  in  Table  (V-l) .   Integration  over  the  wedge  bottom  was 
from  the  apex  to  four  times  the  dump  distance;  this  was  to  insure  that 
contributions  from  the  third  mode  (dumping  about  the  point  x=  3X)  were 
included. 

The  results  of  set  one  are  illustrated  in  Figure  (V-2) .   Note  that 
the  second  maxima  is  severely  depressed  in  the  attenuated  case.   The 
depression  angle  has  been  shifted  slightly  from  11.86  degrees  for  the 
unattenuated  data  to  11.31  degrees  with  Ot  set  to  1.0. 

Figure  V-3  is  a  plot  of  the  results  from  set  two.   All  curves  have 
been  normalized  to  the  initial  slope  for  3  =  1.0.   The  second  and  third 
maxima  decrease  with  increasing  ct,  and  relatively  more  energy  is 
contained  in  the  first  beam.   The  depression  angle  of  the  first  beam 
decreases  slightly  with  increasing  a  .   These  trends  are  summarized  in 
Table  (V-2). 

Figures  (V-4a)  and  (V-4b)  illustrate  the  effects  of  attenuation  on 
different  contributing  modes.   For  Figure  (V-4a) ,  the  source  was  placed 
at  1/3  S.   With  this  source  geometry,  contributions  from  mode  three  are 
suppressed,  and  only  the  first  and  second  mode  contributed  to  the  beam 
pattern  in  the  bottom.   For  Figure  (V-4b) ,  the  source  was  placed  at  ^3 
and  modes  one  and  three  propagate  into  the  bottom  medium.   Both  cases 
were  run  for  a=  1.0,  0.5,  0.1  and  0.0.   As  could  be  expected, 
attenuation  causes  greater  suppression  of  third  mode  interference 
contributions  than  of  second  mode  contributions. 
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Table  (V-l) 
Input  Specifications  for  Attenuation  Runs 


C   =  1500.0  P   =  1.0 

C      =    1578.947  p„   =    1.052632 

2  2 

Cl/C2   =   °"95  Pl/P2   =   °'95 

6    =    5.0°  f   =    150000  hz 

X     =   4.0  Source   distance   =   infinite 

5      =   2.5° 

a    =   0.0,    0.01,    0.1,    0.5,    1.0 


Refer  to  Figure  II-l  for  variable  meaning, 
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Figure  V-2 

Attenuated  and  unattenuated  pressure  amplitudes 
down  from  the  apex 
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Figure   V-3 
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TABLE   7-2 


a 

Y 

Hl 

H2 

Hz/H 

1.0 

12.84 

0.21xl0"4 

NONE 

0.5 

13.39 

.33xl0~2 

.26xl0~2 

1.26 

0.1 

13.32 

0.21 

0.23 

0.90 

0.01 

13.93 

0.54 

0.66 

0.83 

0.0 

13.98 

0.61 

0.74 

0.82 

Y  =  depression  angle  of  first  maxima 


H     =  amplitude  of  first  maxima 
H     =  amplitude  of  second  maxima 
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VI.   ANALYSES  OF  FAR-FIELD  INTERFERENCE  PATTERNS 

The  results  of  the  attenuation  investigation  were  supportive  of  an 
interference  explanation  for  the  multiple  beams.   Subsequently,  two 
analytical  approaches  were  tried  to  predict  the  location  of  interference 
maxima  and  minima  under  unattenuated  far-field  conditions. 

Both  the  Differential  Change  and  Multiple  Lambda  analyses  assume 
very  simple  processes.   Both  assume  point  sources  of  sound  at  the  wedge 
interface  located  at  X  and  3X.      Reallistically  the  sound  sources  behave 
as  phased  arrays  propagating  assymetrically  about  the  dumping  points. 
Phase  differences  between  the  two  dumping  points  are  non-zero.   The 
Differential  Change  formula  becomes  invalid  as  n,  approaches  zero. 
A.   DIFFERENTIAL  CHANGE  APPROACH 

Based  on  the  simple  geometric  considerations  illustrated  in  Fig  (Vl-1), 
the  following  nondimensional  formula  for  the  change  in  interference 
maxima  locations  was  derived. 


An  = 


4    (C  /C  )    sin   9C  sin    $ 


n n    l   (vi-i) 

l[T)2+(l-x)2y/2     fr2+(3-x)2y/2) 
Note  that  equation  (VI-1)  is  dependent  only  on  the  wedge  angle  6  and  the 

sound  speeds  of  the  two  fluid  media. 

A  series  of  WDGBTM  executions  were  made  for  varying  C  /C  and  3,  with 
all  other  variables  held  constant.   Table  (VI-1)  is  a  summary  of  one  case, 
this  case  is  typical  of  results  observed.   Results  will  be  discussed  in 
Section  VI-C. 
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TABLE    VI-1 


Comparison  of   An  from  equation    (VI-1)    and   An  from  WDGBTM  results. 


5.0 


2.5 


1.25 


Cl/C2  =   °-90 

WDGBTM 

An 

An 

maxima 

from 

predicted 

locations 

WDGBTM 

2.38 

1.67 

4.56 

4.05 

1.92 

3.04 

5.97 

2.33 

1.40 

2.33 

3.73 

1.02 

1.60 

4.75 

0.91 

1.38 

5.66 

2.08 

1.22 

1.29 

3.30 

0.66 

0.88 

3.96 

0.49 

0.77 

4.45 
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B.   MULTIPLE  LAMBDA  APPROACH 

The  second  interference  pattern  calculation  also  utilizes  the 
geometry  of  Fig.  (VI-1) .   For  each  C^/C  ratio  and  beta  combination, 
X  and  \2    (the  wavelength  in  the  bottom  media)  was  calculated  according 
to  the  equations 


X  =  ± (VI-2) 

4f  sin   8C  sin   S 

and 

\2   =  C2/f  (VI-3) 


where   is  150khz.   The  ratio  X  /x   =  Q 
is  then  found  and  used  in  the  relation 


NQ   =  \x?+(3-X)2f/2   -  [r)2+(l-X)2]1/2  CVI-4) 

where  W  is  any  positive  integer.   Note  that  use  of  equation  (VI-4) 
assumes  zero  initial  phase. 

To  find  the  expected  maxima  and  minima,  a  plot  of  1  vs  iV  was  made 
for  each  Q.   Maxima  are  expected  for  each  integer  N,   minima  for  each 
N+l/2.   Table  (VI-2)  summarizes  one  case  of  maxima  and  minima  expected 
as  a  result  of  equation  (VI-4)  versus  maxima  and  minima  observed  on 
WDGBTM  output. 
C.   SUMMARY  OF  ANALYSES 

Predictions  from  both  simple  analyses  are  of  the  same  order  of 
magnitude  as  WDGBTM  results.   The  number  of  maxima  and  minima  are 
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TABLE   VI- 2 
OCCURRENCE   OF   MAXIMA   AND  MINIMA 


C1/C2   =0.85 


(DEG) 

T)   MAX   WDGBTM 

n  MAX  EXP 

T)   MIN   WDGBTM 

n  MIN  EXP 

5 

4.03 

2.68 
6.34 

5.39 

4.70 

2.5 

3.59 

2.59 

4.55 

3.76 

5.00 

4.70 

5.66 

5.56 

1.25 

2.02 

2.53 

2.49 

1.76 

2.85 

3.56 

3.16 

3.81 

3.72 

4.77 

4.25 

5.03 

4.68 

5.40 

5.13 

5.62 

5.54 

5.90 

5.94 
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consistantly  predicted  as  increasing  with  decreasing  wedge  angle. 
Occasinal  unexplainable  phenomena,  such  as  increasing  WDGBTM  An  for 
the  3  =  5  on  Table  (VI-1)  and  the  extra  WDGBTM  maxima  for  B  =  5 
Table  (VI-2) ,  were  observed,  but  results  are  generally  supportive  of 
an  interference  effect  occurring. 
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VII.   COMPARISONS 

A.   EXPERIMENTAL  DATA 

Netzorg  (Reference  4)  measured  pressure  amplitude  as  a  function  of 
angle  for  fixed  radii  in  four  media  below  a  wedge  shaped  upper  layer. 
The  angle  domain  version  of  WDGBTM  was  used  to  compare  calculated 
pressure  patterns  with  measured  patterns. 

All  runs  utilized  finite  source  input.   For  each  of  Netzorg 's  four 
experimental  cases,  corresponding  computer  runs  were  made  with  the 
source  distance  set  to  100X,  10X  and  SX,  where  S  corresponds  to  the 
actual  source  distance  used  in  the  tank.   In  all  cases  S  was  less  than 
10X.   The  integrations  were  carried  out  over  one  dumping  distance  and 
four  dumping  distances  for  each  source. 

Figure  (VII-1)  is  typical  of  pressure  amplitudes  in  the  bottom. 
The  source  in  the  experiment  was  only  about    40  wavelengths  from  the 
dumping  point,  even  at  10X  the  source  is  approximately  110  wavelengths 
away.   Since  the  method  of  images  calculations  utilized  to  determine 
the  pressure  amplitude  along  the  wedge  bottom  assume  a  plane  wave,  it 
is  questionable  that  this  program  can  produce  a  realistic  comparison 
with  the  measured  data. 

Figures  (VII-2)  and  (VII-3)  illustrate  the  wide  variation  in  results 
obtained  depending  on  the  source  distance  and  integration  interval . 
Calculations  demonstrate  expected  trends,  i.e.,  pressure  amplitude 
increases  with  decreasing  radius  distance  (CJ)  ,  with  closing  source 
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distance,  and  as  a  function  of  an  increasing  integration  interval  over 
the  wedge  bottom.   The  beam  depression  angle  becomes  more  shallow  as 
source  distance  increases.   Results  at  very  low  angles  are  not 
considered  valid  because  of  assumptions  inherent  in  the  WDGBTM  model. 

The  circled  data  points  of  Figures  (VII-2)  and  (VII-3)  are  a 
transposition  of  data  from  Reference  4.   Since  absolute  amplitude  was 
not  specified  in  the  reference,  rescaling  to  the  case  most  near  to  the 
actual  conditions  was  performed,  i.e.,  the  data  were  normalized  to  the 
case  X     =  l.o  and  srcd  =  3.5X.   The  shape  of  the  experimental  curve  and 
the  apparent  beam  depression  angle  are  the  only  valid  comparisons  that 
should  be  made. 

Comparisons  of  WDGBTM  beam  depression  angles  with  measured  values 
gave  varying  results.   Measured  depression  angles  decreased  with 
increasing  a  ,  WDGBTM  depression  angles  increased.   WDGBTM  results  were 
within  6  percent  to  62  percent  of  measured  values.   No  correlation  was 
observed  between  prediction  accuracy  and  source  distance  or  prediction 
accuracy  and  the  interface  integration  interval.   All  measured  beam 
depression  angles  in  Reference  8   were  for  the  two  shortest  radii  of 
each  case,  the  three  longer  radii  could  only  be  gualitatively  compared 
with  graphical  data.  Possible  attenuation  effects  were  not  considered. 
Numerical  results  are  contained  in  Appendix  D. 

B.   PARABOLIC  EQUATION  APPROACH 

Kuperman  (Reference  9)  describes  a  Parabolic  Equation  model  result 
for  a  wedge  problem,  and  compares  these  results  to  one  set  of  data 
extracted  from  Reference  10. 
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The  Parabolic  Equation  results  show  very  well  defined  beams  in  the 
bottom,  with  no  apparent  interference  effects  as  observed  in  comparable 
WDGBTM  runs  (Figures  V-l  and  V-4) .   Kuperman  used  an  attenuation  factor 
of  0.5dB     per  wave  length,  a  value  equivalent  to  an  input  of  0.984 
Nepers/X  for  WDGBTM.   It  should  be  noted  in  Figure  V-4  that  there  is 
complete  separation  of  the  beams  for  WDGBTM  when  ct  =  1.0. 

Kuperman  states  agreement  within  20  percent  of  a  measured  data  set 
of  Reference  10.   This  measured  data  is  identical  to  Case  1  of  Appendix 
D.   WDGBTM  accuracy  for  this  case  varied  from  zero  percent  error  to  22 
percent  error  with  error  calculated  according  to  the  formula 

,„„    (observed  data)  -  (measured  data) 
100  x  


measured  data 
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VIII.   CONCLUSIONS 

1.  Program  WDGBTM  offers  a  fast,  efficient  means  of  analysing  the 
behavior  of  sound  in  a  fast  bottom  underlying  a  wedge-shaped  medium. 
Effects  of  attenuation  are  easily  demonstrated,  and  modes  desired  are 
readily  controlled  by  specification  of  source  location  and  integration 
interval  selected.   Results  of  the  model  are  not  inconsistant  with 
current  measured  data,  and  illuminate  observations  generated  by  the 
Parabolic  Equation  model. 

2.  A  simple  equation  for  calculating  beam  depression  angle  as  a 
function  of  wedge  characteristics  was  obtained  which  offers  sufficient 
accuracy  for  use  as  a  rapid  analytical  tool.   The  equation  was  derived 

for  acoustic  sources  placed  at  one-half  the  wedge  angle.   Accuracy  degrades 
when  the  wedge  medium  density  to  bottom  medium  density  ratio  decreases 
to  0.5. 

3.  It  is  recommended  that  further  experiments  be  performed  that 
would  further  test  the  model.   New  experiments  should  be  sufficiently 

large  scale  to  ensure  plane  waves  at  the  wedge  bottom  interface  and  to  satis- 
fy  the  k    \r-r    \  >>2ir   criteria.    Bottom  media  that  reasonably  model 
real  ocean  fluid-like  bottom  are  also  recommended. 
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APPENDIX  A 
SAMPLE  INPUT  SET 


o 

o 

o 

o 

sO 

o 

o 

O 

o 

m 


LO  (M 


m 


m 

CVJ 

O 

^> 

CM 

^» 

in 

com 

om 

om 

•  • 

i  • 

»^C\J 

r— i  — * 

o 

o 

• 

• 

i 

i 

O   • 

o  • 

oo 

oo 

oo 

oo 

oo 

ooo 

oo 

0»0 

•in 

•in 

i^-^ 

o-* 

o 
o 


in  in 


CO 

o 

o 

o    o 

f>- 

o 

o 

•      •    m 

J- 

• 

• 

in     w       • 

y-er> 

O 

o 

O 

<*      *-h     ~* 

— »  • 

• 

a.  oo 

1—00 

ZNO 

h- 

»- 

3CMO 

H-      H-      f- 

>—m   • 

3 

3 

0.0    • 

3    r>    ^ 

-*>i- 

CL 

a. 

2-«>J- 

a    a.    q. 

LU 

2 

2 

i— • 

2      2      2 

O 

1— 1 

i-i 

M          M          to* 

QC 

LU 

O       O 

!3 

2 

2 

O 

2    .2    »2 

o 

>— ■ 

*■* 

CC 

'-'O'-'O'" 

00 

< 

< 

O 

<      <      < 

s: 

3: 

o 

21      31      3! 

LU 

a 

o 

oo 

a    o    o 

h- 

o 

o 

a    o    o 

—  O 

LUO 

o 

2   • 

X      i 

=>X      i 

Oh-    • 

•  LU       LU       X       O 

«o     Ot-     i 

OH-      , 

-4—.0      OOJ-iJOt-      o 

s 

a, 

c 
■§ 

0 


§ 

3 


itoooao^aoozooo^omot>aoo 
2in  •v^,lu  •  »lu  •  ••—•in  «^     2  »2  •lu  •  • 

•-■.■HO      QOOQOOt-tO  <-*<OQOO 


44 


H 


p 

Q 

o 

f 

as 

. 

un 

6 

LH 

in 

|-H 

p 

"-H 

M 

Cm 

0 

"* 

Cm 

Cm 

^r 

"3« 

Cm 

< 

lo 

T 

H 

< 

m 

Cm 


< 


CN 


. 

tn 

LO 

Cm 

r»H 

«* 

■*— ' 

s 

< 

*y 

p 
cd 
co 

p 

3 

a 

c 


o 

cd 


cr 
i 

CN 

sr 

c 

0 

c 

•H 

cd 

4-> 

c 

CD 

cu 

0 

3 

0 

•P 
P 

P 
0) 

a) 

Q< 

ja 

u 

rH 

0 

p 

W< 

co 

3 

c 

p 

0 

1 

■H 

G 

w 

(N 

(0 

s 

s 

0) 

UJ^ 

0 

CD 

p 

TJ 

P 

•i-t 

0 

c 

P 

.G 

G 

•r-i 

0 

+J 

•H 

U-l 

U-l 

a 

c 

CD 

U-i 

•H 

^ 

CO 

p 

T3 

0 

CD 

u 

CO. 

c 

P 

P 

3 

0 

CD 

0 

<D 

rH 

V4— 1 

* 

e 

U-l 

1 

n) 

CN 

» 

cd 

> 

Cn 

Q. 

'O 

P 

Di 

3 

G 

0 

3 

C 

•H 

* 

C 

•H 

u 

r-( 

» 

•H 

P 

«. 

p 

a 

4h 

+J 

U^ 

w 

U-l 
0 

Cfl 

rH 

« 

* 

iH 

* 

as 

(N 

CN 

p. 

fll 

CN 

P 

u 

W 

CD 

p 

C" 

0 

,D 

0 

p 

o 

* 

£ 

P 

». 

•H 

rH 

■H 

3 

•H 

H 

Pi 

u 

UJ^ 

a 

J 

er 

G 

0 

•H 

P 

P 

3 

a. 

& 

0 

C 
•P 

CD 

- 

O 

c 

p 

•H 

3 

c0 

0 

£ 

tfl 

0 

CD 

P 

CD 

•H 

rH 

c 

Cn 

•H 

C 

<P 

X 

as 

P 

c 

P 

0 

-H 

O 

IH 

CD 

14-1 

cn 

O 

cn 

c 

G 

G 

•H 

as 

•P 

P 

P 

P 

n 

^£» 

r- 

P 
Cfl 

CN 

m 

T 

Cfl 

■p 

P 
Cfl 

0) 

CD 

0) 

CD 

CD 

CD 

•X3 

C 

C 

G 

rH 

c 

G 

G 

rH 

•H 

•P 

•H 

a) 

•p 

•P 

•P 

CD 

(0 

H 

rH 

rH 

p 

0 

rH 

rH 

rH 

U 
P 

P 
O 

OJ 

CD 

<D 

p 

CD 

CD 

CD 

3 

P 

0) 

CD 

CD 

•p 

CD 

CD 

CD 

0 

■P 

Cfl 

Cfl 

cn 

PI 

Cfl 

Cfl 

Cfl 

Cfl 

PI 

CD 

<3 

CN 

X 

<X> 

LD 
rH 

rH 

in 

<£> 

r- 

„ 

CD 

CD 

CD 

CD 

<D 

rH 

G 

G 

G 

G 

C 

.* 

•P 

•p 

■H 

•H 

•P 

CD 

rH 

rH 

rH 

rH 

rH 

0 

CD 

CD 

CD 

CD 

* 

CD 

CD 

CD 

CD 

CD 

t> 

Cfl 

CO 

Cfl 

Cfl 

Cfl 

CD 

G 
•H 
PI 


in     ko 


co     <n 


^     r-     co 


o 


o 

CN 


<N 
CN 


45 


APPENDIX  B 


Pl/P2   =    '95°° 


Cl/C2 

6 (DEG) 

^DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

15.64 

22.93 

0.68 

6.92 

13.50 

19.08 

0.71 

5.00 

12.02 

17.48 

0.69 

2.00 

9.09 

12.19 

0.74 

1.00 

7.58 

9.54 

0.79 

0.90 

10.00 

19.80 

29.07 

0.68 

6.92 

17.74 

26.66 

0.67 

5.00 

16.33 

21.85 

0.75 

2.00 

11.31 

16.20 

0.70 

1.00 

9.09 

12.73 

0.71 

0.85 

10.00 

24.37 

36.02 

0.68 

6.92 

21.16 

31.05 

0.68 

5.00 

18.42 

26.61 

0.69 

2.00 

13.50 

19.45 

0.69 

1.00 

11.31 

14.36 

0.79 

2L  =  ARCTAN(Z/X)  AT  MAX  AMPLITUDE 

^DEPRESSION 


BANDWIDTH  =  AECTAN(Z/X,     Z/X,  _  ) 


MEAM   BW   =  .73 
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P2/P2=    .9000 


Cl/C2 

B(DEG) 

^DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

15.64 

22.68 

0.69 

6.92 

13.49 

19.08 

0.71 

5.00 

12.02 

17.43 

0.69 

2.00 

9.09 

12.19 

0.74 

1.00 

7.58 

9.54 

0.79 

0.90 

10.00 

19.80 

29.77 

0.66 

6.92 

18.42 

26.75 

0.69 

5.00 

16.33 

21.80 

0.75 

2.00 

12.02 

15.53 

0.77 

1.00 

9.09 

12.90 

0.70 

0.85 

10.00 

24.37 

36.28 

0.67 

6.92 

21.16 

31.22 

0.68 

5.00 

18.42 

26.93 

0.68 

2.00 

13.50 

19.70 

0.68 

1.00 

11.31 

14.95 

0.76 

47 


P2/P2  =    -8000 


Cl/C2 

S(DEG) 

^DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

16.33 

22.39 

.73 

6.92 

14.20 

19.19 

.74 

5.00 

12.79 

17.33 

.74 

2.00 

9.82 

12.13 

.81 

1.00 

7.58 

9.87 

.77 

0.90 

10.00 

20.46 

31.00 

.66 

6.92 

19.14 

26.93 

.71 

5.00 

16.33 

21.90 

.74 

2.00 

12.02 

16.38 

.73 

1.00 

9.82 

12.95 

.76 

0.85 

10.00 

25.03 

36.94 

.68 

6.92 

22.44 

31.72 

.71 

5.00 

19.14 

27.65 

.69 

2.00 

14.20 

20.10 

.71 

1.00 

11.31 

15.32 

.74 
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P2/P2   =    .7000 


Cl/C2 

S  (DEG) 

4- 

DEPRESSION 

BANDWIDTH 

0.95 

10.00 

17.07 

22.24 

6.92 

14.20 

19.34 

5.00 

13.50 

17.17 

2.00 

9.82 

12.13 

1.00 

8.36 

9.65 

0.90 

10.00 

21.16 

32.66 

6.92 

19.80 

26.93 

5.00 

17.07 

22.00 

2.00 

12.79 

16.59 

1.00 

10.59 

13.12 

0.85 

10.00 

25.64 

37.88 

6.92 

23.12 

32.29 

5.00 

19.80 

28.72 

2.00 

14.95 

20.56 

1.00 

12.02 

15.85 

RATIO 


.77 
.73 
.79 
.81 
.87 
.65 
.74 
.78 
.77 
.81 
.68 
.72 
.69 
.73 
.76 
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p      n     =    .5000 


Cl/C2 

8 (DEG) 

4- 

DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

17.74 

22.73 

0.78 

6.92 

15.64 

19.80 

0.79 

5.00 

14.95 

16.91 

0.88 

2.00 

9.82 

12.68 

0.77 

1.00 

9.09 

10.26 

0.89 

0.90 

10.00 

22.44 

35.53 

0.63 

6.92 

21.80 

27.16 

0.80 

5.00 

18.42 

23.70 

0.78 

2.00 

14.20 

17.64 

0.80 

1.00 

11.31 

13.87 

0.82 

0.85 

10.00 

27.20 

42.80 

0.64 

6.92 

25.03 

34.33 

0.73 

5.00 

21.80 

31.00 

0.70 

2.00 

17.07 

21.40 

0.78 

1.00 

13.50 

16.91 

0.80 
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APPENDIX   C 


TABLE  1 

Pl/P2 

C  /C 

r   2 

n 

0.95 

0.95 

0.318 

0.95 

0.90 

0.347 

0.95 

0.85 

0.341 

0.90 

0.95 

0.318 

0.90 

0.90 

0.335 

0.90 

0.85 

0.339 

0.80 

0.95 

0.328 

0.80 

0.90 

0.328 

0.80 

0.85 

0.343 

0.70 

0.95 

0.318 

0.70 

0.90 

0.303 

0.70 

0.85 

0.343 

0.50 

0.95 

0.339 

0.50 

0.90 

0.314 

0.50 

0.85 

0.314 

Average  n 

=   0.329 

Standard  (N-l)  deviation 

=   0.014 

(N-l)  Variance 

=   0.002 
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TABLE  2 

Pl/P2 

0 

Af 

0.95 

10.0 

0.834 

0.95 

6.92 

0.845 

0.95 

5.00 

0.783 

0.95 

2.00 

0.739 

0.95 

1.00 

0.752 

0.90 

10.0 

0.815 

0.90 

6.92 

0.818 

0.90 

5.00 

0.782 

0.90 

2.00 

0.735 

0.90 

1.00 

0.720 

0.80 

10.0 

0.791 

0.80 

6.92 

0.839 

0.80 

5.00 

0.767 

0.80 

2.00 

0.686 

0.80 

1.00 

0.750 

0.70 

10.0 

0.730 

0.70 

6.92 

0.892 

0.70 

5.00 

0.701 

0.70 

2.00 

0.778 

0.70 

1.00 

0.697 

0.50 

10.0 

0.786 

0.50 

6.92 

0.881 

0.50 

5.00 

0.693 

0.50 

2.00 

1.051* 

0.50 

1.00 

0.722 

*Rejected  data 

point 

Average  M 

=  0.772 

Standard  (N-l)  deviation 

=  0.059 

Variance 

=  0.003 
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Table   3 

Cl/C2 

9c(rad) 

0.95 

0.318 

0.95 

0.318 

0.95 

0.318 

0.95 

0.318 

0.95 

0.318 

0.90 

0.451 

0.90 

0.451 

0.90 

0.451 

0.90 

0.451 

0.90 

0.451 

0.85 

0.555 

0.85 

0.555 

0.85 

0.555 

0.85 

0.555 

0.85 

0.555 

6  q 

10.0  -0.217 

6.92  -0.285 

5.00  -0.257 

2.00  -0.215 

1.00  -0.283 

10.0  -0.183 

6.92  -0.232 

5.00  -0.211 

2.00  -0.300 

1.00  -0.317 

10.0  -0.162 

6.92  -0.242 

5.00  -0.255 

2.00  -0.362 

1.00  -0.283 


Average  q     =  -0.254 
standard  (N-l)  deviation    =   0.053 
variation    =   0.0028 
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Table  4 

pl/p2 

Cl/C2 

3 

9BW 

K 

9  DC 

A9D 

0.95 

0.95 

10.0 

15.64 

17.54 

15.35 

-0.29 

0.95 

0.95 

6.92 

13.50 

17.08 

13.60 

0.10 

0.95 

0.95 

5.00 

12.02 

16.94 

12.22 

0.20 

0.95 

0.95 

2.00 

9.09 

17.32 

9.04 

-0.05 

0.95 

0.95 

1.00 

7.58 

18.14 

7.20 

-0.38 

0.95 

0.90 

10.0 

19.80 

16.94 

20.12 

0.32 

0.95 

0.90 

6.92 

17.74 

17.13 

17.83 

0.09 

0.95 

0.90 

5.00 

16.33 

17.55 

16.02 

-0.31 

0.95 

0.90 

2.00 

11.31 

16.43 

11.85 

0.54 

0.95 

0.90 

1.00 

9.09 

16.59 

9.43 

0.32 

0.95 

0.85 

10.0 

24.37 

17.77 

23.61 

-0.76 

0.95 

0.85 

6.92 

21.16 

17.42 

20.92 

-0.24 

0.95 

0.85 

5.00 

18.42 

16.87 

18.80 

0.38 

0.95 

0.85 

2.00 

13.50 

16.72 

13.91 

0.41 

0.95 

0.85 

1.00 

11.31 

17.59 

11.07 

-0.24 

0.90 

0.95 

10.0 

15.64 

17.31 

15.56 

-0.08 

0.90 

0.95 

6.92 

13.49 

16.85 

13.79 

0.30 

0.90 

0.95 

5.00 

12.02 

16.71 

12.39 

0.37 

0.90 

0.95 

2.00 

9.09 

17.08 

9.16 

0.07 

0.90 

0.95 

1.00 

7.58 

17.89 

7.30 

-0.28 

0.90 

0.90 

10.0 

19.80 

16.71 

20.40 

0.60 

0.90 

0.90 

6.92 

18.42 

17.55 

18.08 

-0.34 

0.90 

0.90 

5.00 

16.33 

17.31 

16.24 

-0.09 

0.90 

0.90 

2.00 

12.02 

17.23 

12.02 

0.00 

0.90 

0.90 

1.00 

9.09 

16.36 

9.56 

0.47 

0.90 

0.85 

10.0 

24.37 

17.53 

23.94 

-0.43 

0.90 

0.85 

6.92 

21.16 

17.18 

21.21 

-0.05 

0.90 

0.85 

5.00 

18.42 

16.64 

19.06 

0.64 

0.90 

0.85 

2.00 

13.50 

16.49 

14.10 

0.60 

0.90 

0.85 

1.00 

11.31 

17.35 

11.22 

-0.09 
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Table    4    (continued) 


pl/p2 

Cl/C2 

9DW 

K 

9  DC 

A9D 

0.30 

0.95 

10.0 

16.33 

17.54 

16.03 

-0.30 

0.80 

0.95 

6.92 

14.20 

17.21 

14.20 

0.00 

0.80 

0.95 

5.00 

12.79 

17.25 

12.76 

-0.03 

0.80 

0.95 

2.00 

9.82 

17.91 

9.44 

-0.38 

0.80 

0.95 

1.00 

7.58 

17.36 

7.52 

-0.06 

0.80 

0.90 

10.0 

20.46 

16.76 

21.02 

0.00 

0.80 

0.90 

6.92 

19.14 

17.70 

18.62 

-0.52 

0.80 

0.90 

5.00 

16.33 

16.80 

16.54 

0.21 

0.80 

0.90 

2.00 

12.02 

16.72 

12.38 

0.36 

0.80 

0.90 

1.00 

9.82 

17.16 

9.86 

0.04 

0.80 

0.85 

10.0 

25.03 

17.47 

24.67 

-0.36 

0.80 

0.85 

6.92 

22.44 

17.68 

21.85 

-0.59 

0.80 

0.85 

5.00 

19.14 

16.78 

19.64 

0.50 

0.80 

0.85 

2.00 

14.20 

16.83 

14.53 

0.33 

0.80 

0.85 

1.00 

11.31 

16.84 

11.56 

0.25 

0.70 

0.95 

10.0 

17.07 

17.72 

16.59 

-0.48 

0.70 

0.95 

6.92 

14.20 

16.64 

14.70 

0.50 

0.70 

0.95 

5.00 

13.50 

17.60 

13.21 

-0.29 

0.70 

0.95 

2.00 

9.82 

17.31 

9.77 

0.05 

0.70 

0.95 

1.00 

8.36 

18.51 

7.78 

-0.58 

0.70 

0.90 

10.0 

21.16 

16.75 

21.75 

0.59 

0.70 

0.90 

6.92 

19.80 

17.70 

19.27 

-0.53 

0.70 

0.90 

5.00 

17.07 

16.98 

17.31 

0.24 

0.70 

0.90 

2.00 

12.79 

17.20 

12.81 

0.02 

0.70 

0.90 

1.00 

10.59 

17.88 

10.20 

-0.39 

0.70 

0.85 

10.0 

25.64 

17.30 

25.52 

-0.12 

0.70 

0.85 

6.92 

23.12 

17.61 

22.61 

-0.51 

0.70 

0.85 

5.00 

19.80 

16.78 

20.32 

0.52 

0.70 

0.85 

2.00 

14.95 

17.13 

15.03 

0.08 

0.70 

0.85 

1.00 

12.02 

17.30 

11.96 

0.06 
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Table 

4    (cont: 

Lnued) 

pl/p2 

Cl/C2 

P 

9DW 

K 

9  DC 

MD 

0.50 

0.95 

10.0 

17.74 

16.91 

18.07 

0.33 

0.50 

0.95 

6.92 

15.64 

16.82 

16.01 

0.37 

0.50 

0.95 

5.00 

14.95 

17.90 

14.38 

-0.57 

0.50 

0.95 

2.00 

9.82 

15.89 

10.64 

0.82 

0.50 

0.95 

1.00 

9.09 

18.48 

8.47 

-0.62 

0.50 

0.90 

10.0 

22.44 

16.31 

23.69 

1.25 

0.50 

0.90 

6.92 

21.80 

17.89 

20.99 

-0.81 

0.50 

0.90 

5.00 

18.42 

16.82 

18.86 

0.44 

0.50 

0.90 

2.00 

14.20 

17.53 

13.95 

-0.25 

0.50 

0.90 

1.00 

11.31 

17.54 

11.10 

-0.21 

0.50 

0.85 

10.0 

27.20 

16.85 

27.80 

0.60 

0.50 

0.85 

6.92 

25.03 

17.50 

24.62 

-0.41 

0.50 

0.85 

5.00 

21.80 

16.96 

22.13 

0.33 

0.50 

0.85 

2.00 

17.07 

17.96 

16.37 

-0.70 

0.50 

0.85 

1.00 

13.50 

17.84 

13.03 

-0.42 

Average  K  =  17 .22 

Standard  (N-l)  deviation   =   0.51 

(N-l)  variance   =   0.26 


56 


Case    1 
2.12 


10 


100 


APPENDIX  D 

X 

Measured 

Measured 

max 

a 

9DW 

BW 
w 

9D 

BW 

1 

1.60 

24.5 

10.33 

1 

1.10 

22.0 

11.34 

1 

0.8  3 

19.0 

12.88 

1 

0.66 

16.5 

15.27 

13.5 

15 

1 

0.55 

14.5 

18.38 

15.5 

15 

4 

1.60 

23.0 

6.28 

4 

1.10 

22.1 

8.19 

4 

0.83 

19.0 

12.20 

4 

0.66 

16.0 

13.88 

13.5 

15 

4 

0.55 

15.0 

14.71 

15.5 

15 

1 

1.60 

17.5 

19.26 

1 

1.10 

16.0 

20.29 

1 

0.83 

16.5 

20.90 

1 

0.66 

13.5 

22.41 

13.5 

15 

1 

0.55 

12.5 

23.19 

15.5 

15 

4 

1.60 

16.0 

13.85 

4 

1.10 

13.5 

14.96 

4 

0.83 

13.0 

15.36 

4 

0.66 

12.0 

15.67 

13.5 

15 

4 

0.55 

12.0 

16.33 

15.5 

15 

1 

1.60 

16.5 

22.68 

1 

1.10 

15.5 

23.27 

1 

0.83 

14.0 

23.91 

1 

0.66 

13.0 

24.81 

13.5 

15 

1 

0.55 

11.5 

26.00 

15.5 

15 

4 

1.60 

12.0 

15.18 

4 

1.10 

12.0 

14.93 

4 

0.83 

11.5 

15.25 

4 

0.66 

11.0 

15.94 

13.5 

15 

4 

0.55 

10.5 

16.95 

15.5 

15 

Cl/C2  = 

0.8854 

P-./P 

2  =  0.8415       6 

=  1.52° 
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Case  2 

X 

s 

X 

max 

a 

^W 

BW 
w 

Measured 
9D 

Measured 
BW 

3.50 

1 

1.21 

21.5 

18.31 

1 

0.83 

17.5 

20.60 

1 

0.63 

14.5 

24.19 

1 

0.50 

12.5 

26.62 

13.3 

15.5 

1 

0.42 

10.0 

27.41 

14.0 

14.33 

4 

1.21 

22.5 

15.04 

4 

0.83 

16.0 

19.24 

4 
4 
4 

0.63 
0.50 
0.42 

14.0 
13.0 
11.5 

19.  32 
20.34 
22.75 

13.3 
14.0 

15.5 
14.38 

10 

1 

1.21 

18.0 

24.19 

1 

0.83 

14.5 

25.13 

1 
1 
1 

0.63 
0.50 
0.42 

14.5 

12.0 

9.5 

28.69 
29.56 
30.18 

13.3 
14.0 

15.5 
14.38 

4 

1.21 

16.0 

16.07 

4 

0.83 

14.0 

17.97 

4 
4 
4 

0.63 
0.  50 
0.42 

14.0 
12.5 
11.5 

20.43 
22.70 
24.96 

13.3 
14.0 

15.5 
14.38 

100 

1 

1.21 

18.5 

26.64 

1 

0.83 

16.0 

27.72 

1 
1 
1 

0.63 
0.50 
0.42 

14.0 

11.5 

9.5 

29.29 
30.89 
31.55 

13.3 
14.0 

15.5 
14.38 

4 

1.21 

14.0 

17.65 

4 

0.83 

13.5 

17.94 

4 
4 
4 

0.63 
0.50 
0.42 

12.5 
11.5 
11.0 

19.37 
21.69 
23.80 

13.3 
14.0 

15.5 
14.38 

C. 

L/c2  = 

0.8854 

V 

3   =  0.8415       6 
2 

=  2.51° 
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Case  3 

9DW 

BW 
w 

Measured 

Measured 

s       x 

max 

a 

9D 

BW 

6.15 

1 

0.58 

13.5 

30.13 

1 

0.40 

8.0 

32.50 

1 

0.30 

7.0 

33.47 

1 

0.24 

7.0 

35.75 

16.4 

25.5 

1 

0.20 

7.5 

38.56 

20 

21 

4 

0.58 

14.0 

24.62 

4 

0.40 

10.5 

30.65 

4 

0.30 

10.5 

35.06 

4 

0.24 

10.5 

41.51 

16.4 

25.5 

4 

0.20 

11.0 

20 

21 

10 

1 

0.58 

14.5 

24.60 

1 

0.40 

12.5 

30.47 

1 

0.30 

11.5 

35.89 

1 

0.24 

11.5 

41.62 

16.4 

25.5 

1 

0.20 

12.0 

20 

21 

4 

0.58 

14.5 

33.94 

4 

0.40 

9.0 

35.10 

4 

0.30 

7.5 

36.04 

4 

0.24 

7.5 

38.04 

16.4 

25.5 

4 

0.20 

7.5 

40.78 

21 

21 

100 

1 

0.58 

14.0 

35.19 

1 

0.40 

8.5 

35.83 

1 

0.30 

7.0 

36.47 

1 

0.24 

7.0 

38.37 

16.4 

25.5 

1 

0.20 

7.5 

40.82 

20 

21 

4 

0.58 

13.5 

25.14 

4 

0.40 

11.0 

30.32 

4 

0.30 

10.0 

35.38 

4 

0.24 

10.0 

40.97 

16.4 

25.5 

4 

0.20 

11.0 

35.02 

20 

21 

Cl/C2 

=  0. 

8854 

P1/p 

2  =  0.8415       6 

=  4.41° 
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Case    4 


*s 

*max 

a 

e 

DW 

BW 
w 

Measured 
9D 

Measured 
BW 

3.10 

1 

1.31 

18.5 

15.88 

1 

0.90 

14.5 

17.62 

1 

0.68 

13.0 

20.75 

1 

0.54 

10.0 

22.84 

12.67 

12.88 

1 

0.45 

8.0 

24.58 

14.7 

13.5 

4 

1.31 

17.5 

9.65 

4 

0.90 

15.5 

13.38 

4 

0.68 

14.0 

17.32 

"4- 

0.54 

12.0 

21.72 

12.67 

12.88 

4 

0.45 

11.0 

25.21 

14.7 

13.5 

10 

1 

1.21 

17.0 

21.32 

1 

0.90 

14.0 

22.55 

1 

0.68 

12.0 

24.35 

1 

0.54 

10.0 

26.14 

12.67 

12.88 

1 

0.45 

8.0 

26.62 

14.7 

13.5 

4 

1.31 

14.5 

13.80 

4 

0.90 

13.0 

15.71 

4 

0.68 

12.0 

17.96 

4 

0.54 

11.5 

20.56 

12.57 

12.88 

4 

0.45 

10.0 

22.32 

14.7 

13.5 

100 

1 

1.31 

16.0 

22.58 

1 

0.90 

13.5 

23.71 

1 

0.68 

11.5 

25.18 

1 

0.54 

9.5 

26.56 

12.67 

12.88 

1 

0.45 

7.5 

26.76 

14.7 

13.5 

4 

1.31 

12.5 

15.05 

4 

0.90 

11.5 

15.87 

4 

0.68 

11.0 

17.27 

4 

0.54 

10.0 

19.32 

12.67 

12  .88 

4 

0.45 

9.0 

21.15 

14.7 

13.5 

C^/Cj    =    0. 

,9214 

P1/  P 

■2    =    0.9  07  3                  3    =    2.6  6° 
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